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Abstract. The one-loop HQET heavy-heavy vertex diagram with arbitrary powers of all three denominators and 
arbitrary residual energies is investigated. Various particular cases in which the result becomes simpler are considered. 



1 Introduction 

The heavy-heavy quark current in HQET (see, e.g., 0~|[2]]) trans- 
forms a heavy quark with velocity v\ into a heavy quark with 
velocity w 2 . Loop diagrams with a velocity-changing vertex ap- 
pear, e.g., when the anomalous dimension of the heavy-heavy 
current O or correlators involving this current (see [1J) are 
considered. They can be reduced to a set of master integrals. 
Some of them have the form of the one-loop vertex diagram 
with various e-dependent powers of the denominators. 

Here we consider this one-loop diagram (Fig. [T} in dimen- 
sional regularization in the most general case. It is 
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It depends on the residual energies cji,2 = Pi, 2 ■ vi,2 and the 
Minkowski angle coshi? = v\ ■ v 2 . It is symmetric: 



I(n ll n 2 ,n 3 ;-d;uj ll uj2) = J(n 2 , n x , n 3 ; cj 2 , w x ) 



(2) 



in the following, we shall not explicitly write down relations 
obtained by this obvious symmetry. This diagram vanishes at 
integer n 3 < 0. For integer n 2 < 0, it becomes the self-energy 
diagram with a numerator; if, in addition, n\ is integer and 
ri\ < 0, the diagram vanishes. The diagram has cuts in 0*1,2 
from to +00. We shall consider 0^2 < 0, other cases can be 
treated by analytical continuation. 

In physical applications, 711,2 = mi,2 + %h,2£> n 3 — m 3 + 
Z 3 e, where m t , U are integer, and U > 0. All integrals with 
given li can be reduced, using integration by parts (4), to 3 
master integrals 



1(1 + 2he, 1 + 2l 2 e, 1 + he; tf; wi, w 2 ) 
I(2he, 1 + 2l 2 s, 1 + he; d; uj u uj 2 ) , 
1(1 + 2he, 2l 2 e, 1 + Z 3 e;tf;wi,o;2) . 



(3) 
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Fig. 1. The one-loop HQET heavy-heavy vertex diagram in momen- 
tum and coordinate space 



If l\ — 0, the second integral is trivial: 

1(0, n a , n 3 ; J?; ui,uz) = I(n 2 ,n 3 )(-2u 2 ) d - n2 ~ 2n3 , 
r(m + 2n 2 - d)r(d/2 - n 2 ) 
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see (similarly for the third integral if l 2 — 0). Reduction 
algorithms for all cases has been constructed and implemented 
in REDUCE, and can be downloaded from [6\. 

Using HQET Feynman parametrization (see, e.g., J3J), we 
can write the diagram ([l} as an integral in 2 parameters j/12 
which have dimensionality of energy and vary from to 00. 
After the substitution y\ = yx, y 2 = y(l — x), the integral in 
y can be calculated: 



I(ni,n 2 ,n 3 ;-d;uj 1 ,Lo 2 ) = /(ni +n 2 ,n 3 ) 
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(5) 



This integral depends on d and n 3 only via I; therefore, shifting 
d by ±2 (7) is equivalent to shifting n 3 by ±1. It is Lauricella 
function F D (8) 

7(ni,n2,n 3 ;i?;a;i,a; 2 ) = I(nx + n 2 , n 3 )(~2uj 2 )~ r 
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x F D (n 1 ;l,l,r;n 1 +n 2 ;l- e^ 7 l- e (6) 



where 



UJ 2 



Substituting the explicit form of Gegenbauer polynomials 

[fc/2] 



(7) 



( l) m (0fc-m (2t) k ~ 2m 
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This result can also be obtained in coordinate space (Fig. [TJ 
Substituting t\ = tx, t 2 = t(l — x), we can calculate the inte- and interchanging the order of summations, we get 
gral in t and reproduce (0. Another one-dimensional integral 
representation can be obtained by separating the k space into 
the 2-dimensional longitudinal and (d — 2) -dimensional trans- 
verse subspaces [9|, but this representation seems more com- 
plicated! 

The case v\ — v 2 = 0) has been considered in iPPJll : 



= E E 
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Now we substitute k = p + m and interchange the order of 
summations again: 



I{n 1 ,n 2 ,n 3 ;0]uj 1 ,uj 2 ) = I{n\ + n 2 , n 3 )(-2u 2 ) r 
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If ni 2 are integer, it reduces to trivial cases (0]). 

In the single-scale case cji = cj 2 , the integral (f5J) gives the 
Appell function Fi M 



p=0 r ^ m! (P- m ) ! ( n l+^2) m +p+ ;i 

Here the inner sum is a terminating hypergeometric series 
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I(ni,n2,n 3 ;'d;u),u}) = I(ni + n 2 ,n 3 )(-2w) r 
x Fi Oi; Z, J; ni + n 2 ; 1 - e*, 1 - e^) . (9) expressing it via P functions, we have 



Also, 



I(ni,ra 2 ,7i 3 ;i?;a;,0) = I(d-n 2 -2n 3 ,n 2 ,n 3 ;r);uj,uj) (10) 
(similarly for lo\ = 0); so, these cases reduce to oj-y = lu 2 . 

2 Exact results 

In the Euclidean region # = i-&E, let's consider 

7(ni,n 2 ,n 3 ;i^;wi,w 2 ) 
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Alternatively, interchanging the order of summations in I (fT~4T >. 
we can write it as 



Let's denote 
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A = 1 - 2tz + 2 , i = sin— , z = 2ta. (12) 



Unfortunately, these sums seem not to be expressible via any 
well-known special functions. 

The result simplifies if m = n 2 . Using the quadratic trans- 
formation 



If t < 1/2, we can expand 

00 
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where C l k (t) are Gegenbauer polynomials; results for t > 1/2 
can be obtained by analytical continuation. Expanding also J?/w 2 
in x, we obtain from (|5]l 



and expanding the new hypergeometric function, we can ex- 
press / via the Appell function F 3 |8|: 
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Here (Y + l)/2 + Z = 711 + 1/2, and we can reduce this function 
to the Appell F\ using 



1 In the ordinary relativistic theory, the one-loop vertex diagram in 
the most general case is given by Fd 1 1 1 ; the massless one via Appell 
F4 functions 1111 : and the one with any masses but all powers of the 
denominators equal to 1 via Appell F\ 1121 . 
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Finally, we arrive at the result 



J(ni,ni,7i3;i9;wi,£J 2 ) = 7(2ni,n 3 ) (-2^Mcj 2 ) 
„ (, r 1 (wi-w 2 ) 2 1-coshtA 

At -d = it reduces to (0, due to the quadratic transformation 
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Another case when the result simplifies is u)\ = u> 2 . Then 
the only non-vanishing term in the sum (TBl i is s = 0, and 
from ( fT3T l we obtain the result 



I(ni,n 2 ,n 3 ;d; = 7(m + n 2 , n 3 )(-2u) r 

( ni,n 2 , 1 1 - coshiT 

x 3^2 I m+n 2 m+112 + l 
V 2 ' 2 
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This form is much simpler than (|9j. Due to ( TTOb . we also know 
I(n\, n 2 , n 3 ; cj, 0). If n\ — n 2 , this hypergeometric function 
reduces to 2 Fl, which agrees with (I17t . 

A simple result for ni = n 2 = 1 can be derived directly 
from (0. Substituting x = [-1 + (e" 5 + l)z] /(e^-l), we get 



/ = (e 2 " - 1) 
F ' 6 
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where 



f z-\l-z)- l {l-yz) 2l - 2 dz, (19) 
Jo 
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Expanding the two brackets in ( [T9| > and integrating, we obtain 
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This Appell function can be reduced using 

Ft(a;&i,& 2 ;&i + &2;^i,2 2 ) = 
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and so 
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At last, we arrive at 
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At wi = lo 2 this should be equivalent to (jTSj. 

3 Expansions in e 

The hypergeometric wi = u) 2 result ( TT8l can be expanded in 
e. Expansion of such hypergeometric functions was considered 
in a number of papers [ 14 1. It can be performed automatically 
using the Mathematica package HypExp lfT51 . Coefficients 
are expressed via harmonic polylogarithms 1161 . and can be 
simplified by the package HPL ifTTl . It is better to work in Eu- 
clidean space, where arguments of logarithms and polyloga- 
rithms are away from branch cuts, and to do analytical contin- 
uation to Minkowski space at the end. 

All integrals reduce to 3 master ones. Their expansions have 
the form 

7(1 + 2he, 1 + 2l 2 e, 1 + Z 3 e; u, u) 
1(2 + 2{h + l 2 )e, 1 + l 3 e)(-2u)-^+ l ^+ 1 ^ 

= ^^ + a ^ £ + a ^ £2 + --^ > W 

7(1 + 2he, 2l 2 e, 1 + Z 3 e; to, lo) 

7(2 + 2{h + Z 2 )e, 1 + l s e)(-2u) 1 -^ 1 +^+ l '+ l ») E 

= 1 - 2l 2 e [r& + h (0)e + b 2 (i?)e 2 + ■ • •] , (22) 

where r = tanh(i9/2) (expansion of 7(2Zi£, 1 + 2l 2 e, 1 + 
Z 3 e;$;u;,u;) is given by the formula symmetric to (l22l ). At 
l 2 = 0, the right-hand side of (l22l is exactly 1; a n ($) are odd 
functions (~ i? 3 at $ — > 0); 6„(i?) are even functions (~ i3 2 at 
-> 0). 

The first corrections are 



b 1 {d)=TC 1 {d)-\j'd 2 , 



(23) 
(24) 



ci = —2l\ 2 L 2a — jL 2 b , 
where l 12 = h + l 2 , j 3 = 1 + h, 3 = h - h%, f = 33 ~h + h, 
L 2a = Li 2 (r) - Li 2 (-r) , 
L 2b = Li 2 I - Li 2 ( I Li) . 



L = - 2 g 4 = lo § [ 2 cosh 2 
Two dilogarithms are not independent: 



Li 2 



1 +T 



Li 2 



1 - T 



= -7^ + — + — : 
4 6 



7 2 f, is written here in the manifestly odd form. 
The second corrections are 



02 W =c 2 {#)-4l 2 2 L 2a -2l 12 jL 



2b ■ 



(25) 



b 2 (ti)=TC 2 (#)+4j'j 3 L 3d 

- 2 [(3/i + Z 2 )j3 - Z 2 + Z 2 ] 7 3e 

- h-0 [4/i 7 2a + (j 3 - 3/i + Z 2 )7 2b + /7tf] , (26) 
c 2 (i?) = 4Z i2 j 3 7 3a + 2jj 3 L 3b - 2l 12 jL 3c 
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where 

L3a = Li 3 (r) - Li 3 (-r) , 

^(p)-MrS). 

L 3d = Li s (^)+Li 3 (j-fj 

--L 3 + ^L-I C (3), 
3 6 4^ v 

L -- Li3 (iT7) +Li3 (r^) 

{Lzd can be written via a single trilogarithm Li 3 (— (l±r)/(l=F 
r)). 

The next correction can be derived by HypExp; it contains 
harmonic polylogarithms which cannot be reduced to ordinary 
polylogarithms. 

If rii = n-x = 1, I ( fTTT i depends on n 3 and e only via 
the product j 3 e, see ©. The coefficients ( l23b . ( fZSb satisfy this 
requirement, and coincide with the expansion of < f20b at u)i = 

UJ 2 . 

At small d, keeping only the leading d 3, terms in 0^2(1?) 
and the leading d 2 terms in bi t2 ("&), we reproduce the d 2 term 
in ( TT~8b expanded in e. 



4 Conclusion 

In the single-scale case uj\ — uj 2 , the vertex diagram with 
arbitrary powers of all denominators is given by the surpris- 
ingly simple formula ( fT~8b . The cases wi 2 = reduce to this 
one ( fTOb . Expansions in s has been derived in Sect. [3] If lo\ ^ 
u>2, we were able to obtain relatively simple formulas for n\ = 
n 2 ( E) and especially n\ = 112 = 1 d20b . 

The integral (Q3 with m = ri2 = n 3 = 1 was considered 
in (T8 | for general wi^ (in fact, even a more general integral 
with a nonzero mass m of the light line was calculated up to 
O(e )). Our 01 (|23]l with Zi = l 2 = I3 — reproduces the 
formula (37) from this paper. 

Integrals with m = n% = n 3 = 1; n± = n 2 = 1, n 3 = 
1 + e; ni = 1 + 2e, n 2 = n 3 = 1 were considered in the 
Appendix of OH (the formulas (69), (71), (72)). It is difficult 
to make sense of these formulas, because all integrals in the 
left-hand sides of equations in this Appendix are equal to 0, 
and the right-hand sides have dimensionality different from the 
left-hand sides. Evidently, the author assumed some different 
integrals in the left-hand sides, and some dimensional factors in 
the right-hand sides, but it seems impossible to guess their real 
meaning. In particular, these formulas cannot be interpreted as 
our integrals with u)x = uj 2 with appropriate powers of — 2u 
inserted into the right-hand sides, because already their val- 
ues at d — differ. It would be interesting to compare the 
formulas (69), (71) (containing 2 Fi hypergeometric functions) 
and (72) (containing the Appell F\ function) to our results. The 
O(e ) term in (69) contains dilogarithms similar to those in 



ai (|23]l (at h = l 2 = h = 0), but the leading C(e _1 ) terms 
differ. 
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